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Abstract. Stochastic cooling of trapped atoms is considered for a laser-beam configuration 
with beam waists equal or smaller than the extent of the atomic cloud. It is shown, that various 
effects appear due to this transverse confinement, among them heating of transverse kinetic 
energy. Analytical results of the cooling in dependence on size and location of the laser beam 
are presented for the case of a non-degenerate vapour. 
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1. Introduction 



Cooling techniques for atoms play a crucial role in modern physics. They have allowed 
for the localisation of atomic gases by weak trapping forces and step by step have enabled 
the reach into the domain of ultracold gases. At such low temperatures the peculiar 
quantum-statistical properties of the atoms become immanent in various effects of cold 
atomic collisions. Bose-Einstein condensation [1-5] and the recent production of Fermi- 
Dirac degenerate gases [6-9] are limiting cases of the now existing experimental feasibilities. 
Moreover, the implementation of atom-lasers [10-12] and microstructured traps on so called 
atom chips [13-16] show the vast potential of applications. 

The typical strategy to generate a Bose-Einstein condensate from a moderately cold 
sample of trapped atoms is to apply different cooling techniques in sequence: it usually starts 
with laser cooling [17-19] and ends with evaporative cooling [20-22]. The latter technique 
seems to be unbeaten as of yet for the final cooling step. Laser cooling, that relies on cycling 
transitions where photons are spontaneously emitted, does not provide the ultimate cooling 
power, due to the reabsorption and scattering of the emitted photons. However, the drawback 
of evaporative cooling is well known to be its intrinsic loss of atoms. Since hot atoms are 
released from the trapping potential to reach a colder sample, a substantial atom loss has to 
be taken into account that ultimately limits the size of the condensed sample. Furthermore, 
as does sympathetic cooling [23], it requires sufficiently strong atomic collisions for thermal 
re-equilibration. 

Given a prepared sample of condensed atoms, a multitude of technical and possibly 
fundamental noise effects lead to a finite lifetime of the condensate state. Among these noise 
effects there are collisions with background vapour, electromagnetic noise sources via the 
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trapping potential, scattering of light, etc. The study of these detrimental disturbances and the 
development of methods to reduce their impact on the condensate will be a challenging task 
for the future. One way to compensate for such heating effects may be simply the continuous 
application of cooling during the entire experiment. This may partially compensate the 
heating and thus extends the lifetime of the condensate. The only technique working so far 
at these temperatures is evaporative cooling. Its continuous application, however, would be 
rather unfortunate for the condensate, since though the lifetime of the condensate may be 
extended, its size in terms of atom number will continuously decrease. 

Some years ago, Raizen et al. have proposed the use of stochastic cooling for trapped 
atoms [24]. It is a successful method in high-energy physics [25,26] where the transverse 
motion of a particle beam has to be collimated and cooled. Clearly, the energies involved 
there are not in the regime important for an application to trapped atoms. However, classical 
numerical simulations have shown the feasibility of stochastic cooling also for trapped atoms 
[24,27]. Furthermore, it has been recently shown, that also at ultralow temperatures this 
technique reveals cooling [28]. Thus it may perhaps be utilized to stabilise an atomic Bose- 
Einstein condensate. 

It should be pointed out that on the single-atom level feedback control of atomic position 
has been theoretically studied [29, 30] and experimentally realised in optical lattices [31] and 
high-quality cavity fields [32]. 

In this paper we extend our analysis of stochastic cooling of trapped atoms to include 
also effects due to the transverse confinement of atoms. The latter has its origin in the finite 
beam waist of the employed control-laser beam. At temperatures above the condensation 
point we give analytic results of the cooling and discuss its optimisation with respect to size 
and location of the control-laser beam waist. 

In Sec.|2]we explain the method of stochastic cooling of atoms and derive the expression 
for the single-atom density matrix after the single cooling step. Given this result at hand, 
in Sec. [3] we calculate the total energy change of atoms due to a single step of stochastic 
cooling in terms of quantum-statistical averages. In Sec. |4] the regime of a non-degenerate 
gas is considered for which analytical expressions for the energy contributions are obtained. 
Moreover, the dependence of cooling on geometrical parameters is discussed. Finally, in 
Sec.[5]conclusions are given. 

2. Stochastic cooling 

The method of stochastic cooling of trapped atoms consists of the repeated application of two 
operations: the measurement of the momentum of atoms and the subsequent application of a 
kick to compensate for the measured momentum. Several aspects are important and should 
be emphasised for an understanding of the working of this technique. First of all it is not done 
on a single atom but on a large set of atoms. Those atoms that are subject to measurement 
and kick are specified by their spatial location in a given volume of space. In the experiment 
that volume is defined by the spatial extent of the laser beams that implement the required 
operations. 

Since in the experiment, at the time of measurement of momentum, it is usually unknown 
how many atoms contributed to the measured signal, the momentum per atom averaged over 
the atomic ensemble is not the measured observable. To obtain this momentum per atom 
one would in fact need knowledge on the precise number of atoms, that contributed to the 
measured signal. What can, instead, be assessed by the measurement is the total momentum 
of the atoms. This is the sum over the atomic momenta, since each atom equally contributes 
to the signal. 
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Given the measured total momentum of the set of atoms, for compensating it, a (optical) 
field is turned on to provide the necessary kick by its interaction with the atoms. Since each 
atom separately interacts with the field, one can only apply a common kick to each atom. 
The determination of the required kick per atom necessarily involves a characterisation of the 
number of atoms in the set, given that only the total momentum is known. Since the atom 
number will not be measured, a priori information is required for estimating the actual atom 
number. Clearly, this way atom-number fluctuations, whether classical or quantum in nature, 
cannot be coped with, which shows an intrinsic source of imperfection of the method. 

Furthermore, since a measurement on a single system and not a series of measurements 
on identically prepared systems is performed, not ensemble averages are measured. 
Depending on the measurement resolution strong correlations between the atoms are induced 
by the measurement projection, since a huge number of microstates of the atoms may 
be associated to the same observed measurement outcome, which form a complicated 
superposition state. 



2.1. Single-atom density matrix 

Here we consider a full three-dimensional model and consider spatial confinement of the 
volume where atoms are manipulated, see Fig.[0 For simplicity we assume only one laser- 
beam profile, despite the fact that several laser beams are involved in the implementation 
of the required operations [33-35]. The control-laser beam is directed along the z-axis 
and its transverse profile in x and y directions is described by the beam-waist function 
w±_{v) = w±(x, y). Thus, the z component of the total momentum of atoms inside the beam 
P w is measured and then compensated to zero by means of a negative feedback loop. Using 
the atomic field operator </>(r) for bosonic atoms, i.e. with commutator 

[0(r),0t( r ')] =( 5( r _ r ') ) (1) 
this observable can be written as (H= 

P w =-i [ dVw X {T)ft(T)dJ(T). (2) 



The many-body quantum state of the atomic cloud after a single operation of stochastic 
cooling can be given as an integral over all possible measurement outcomes P for the 
measured total momentum: 



Q- 



= J dPU{P)M{P)Q^M\P)U\P). (3) 

In this expression the initial many-body density operator is denoted as Q- and the final 
one, after the single feedback operation, is denoted as q + . It is assumed here, that the 
measurement and shift of momentum can be performed on a time scale At much faster 
than the characteristic dynamics of the free system, i.e. At <C u^ 1 with u> being the trap 
frequency. Then measurement and shift can be taken as instantaneous processes without time 
delay between them. 

The measurement of momentum P w with specific outcome P is described by the 
resolution amplitude [36, 37] 



J) Throughout the paper we use H = 1. 
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Figure 1. Geometry of the feedback setup. The laser beam is aligned along the z axis, it 
determines the size of the feedback region. 



where a denotes the measurement resolution. Applied on a momentum eigenstate \Pq) it gives 
the probability amplitude to observe the value P.f f A quasi canonically conjugate operator 
for the centre of mass of the atoms in the laser beam, being experimentally accessible, can be 
defined analogously as [39] 

Q™ = jj- fdVw L {r)ft{v)z^v). (5) 

Since N e is an estimated atom number, the commutator relation of operator P w and Q w 
reveals a deviation from the usual canonical form: 

[Q w ,P w ]=iN w /N e . (6) 

Here the true atom-number is defined as the operator 

N w = ( dVwl{Y)ft{v)4>{v). (7) 



The modified commutation relation has an impact on the action of the shift operator U(P), 
that is supposed to produce a shift of P w by — P, and that is defined as 

U(P) =exp(-iPQ w ). (8) 
Transforming the observable P w by use of to the unitary transformation (|8} we obtain 

U\P) P w U{P) =P W - PN W /N e . (9) 
Thus an optimal shift by — P is produced only on average when additionally estimating 

N e = (N w ). (10) 

Nevertheless, atom-number fluctuations will always deteriorate the perfection of the shift 
operation. In view of the lack of knowledge on the true atom number N w , the estimate dlOt 
represents an optimum. Thus in the following we use this justifiable estimate. 

ff The set of operators Aft (P)M(P) forms a positive operator-valued measure [38]. 
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For our further derivations it is convenient to calculate the single-atom density matrix 
from Eq. l|3}- It is defined as 

t 7(r 1) r 2 ) = (^(r 2 )0(n)) J (11) 

and can be calculated as a trace over the many-body density operator given in Eq. 0. In this 
way the single-atom density matrix after (+) an operation of stochastic cooling reads 

o+(n,r a ) = jdP^{P)U\P)^{v 2 )^{v l )U{P)M{P)^_, (12) 

where (...)- denotes tracing over the many-body density operator that represents the 
quantum state before the feedback operation. 

The action of U (P) on a field operator results as a c-number exponential factor 

fr(P)4>(r)U(P) = fa) exp 

Moreover, using the Fourier representation of the resolution amplitude 

M(P) = f dqM(q) e ^ p - p '\ 



-izw ± {r)P/{N w ) 



with 



M(g)=\/^exp(-aV) 
the single-atom density matrix can be rewritten as 

o + fa.ra) = [dp[dq[dq'( e** p <° & ^fai) e~*' p « 

Z 2 W_l(v 2 )-ZiW1_{yi) , 



M*(q)M(q') 



x exp< iP 



This result can be further simplified using the transformation 



exp 



iqPu 



(r) exp 



-iqP u 



which results in 

tr+(ri ) r 2 )= [ dP I dq I dq' M*{q)M(q') ex-pi iP 



x,y, z-qwx(r)), 

z 2 w± (r 2 )-ziwj_(ri) 

(N w ) 

x (ft (x 2 ,y 2 , z 2 -qw±(r 2 )) 4>(xt, y x , Zx-qw^Tx)) e l(q ^< 
Performing then the P and q 1 integrations we finally obtain 

a+(r 1; r 2 ) = 2n J dqM*{q)M(q+[z 1 w±{r 1 )-z 2 w ± {r 2 )]/(N w )) 

x (ft(x 2 ,y 2 ,z 2 -qw ± (r 2 ))$(x l7 y 1 ,z 1 ~qw ± (r 1 )) ^N^)-^^)]/^,) 



(13) 



(14) 



(15) 



(16) 



(17) 



(18) 



,(19) 

which shows that the single-atom density matrix after the feedback depends in general on 
higher atom-atom correlations before the feedback. In Eq. dl9l this is encoded by the 
occurrence of the exponential operator. 
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3. Feedback-induced change of energy 

Important features of the application of a single feedback step can be given in analytical form. 
For example, the difference of energy after and before the feedback step AE can be calculated 
from Eq. il9i . The information contained therein allows us to recognise noise sources and 
determine optimal parameters for maximum cooling. The parameters that can be optimised 
are the measurement resolution a and the geometrical characteristics given by the size and 
location of the control-laser beam with respect to the trapping potential. 

In the following we use the thermal equilibrium state of the atomic ensemble to calculate 
the average energy change. This allows us to obtain a natural description of cooling in terms 
of energy AE(T) that is subtracted or, possibly, added by a feedback operation at a certain 
temperature point. However, this approach does not necessarily reflect the most general 
experimental situation, since specific correlations generated step by step in the feedback 
process are not taken into account. These correlations may possibly lead to enhancement 
of cooling via various effects, as shown in Refs [24,27]. Thus the results of this paper, where 
we assume a thermal equilibrium state, represent the leading cooling/heating mechanisms for 
the quasi-equilibrium case. 

Having the expression for the single-atom density matrix fl!9l > the change of total energy 
due to the application of a single feedback step can be formulated. The Hamiltonian of the 
system of non-interacting atoms in the isotropic, harmonic trap potential is 

' ] kr), (20) 



H = JdV ft(r) 



- H r 2 

2m 2 



where m and u> are the atomic mass and the vibrational trap frequency, respectively. We 
divide this Hamiltonian into parts describing the energy of the motion in z direction, i.e. in 
longitudinal direction with respect to the measured momentum P w , and into parts related to 
the energy of the transverse motion in the xy plane. 

The average total energy of the system in a given many-body quantum state can be 
written as the sum of these different contributions as 

E=(H) = T,| + T± + Vj, + V±. (21) 

Using the definition of the single-atom density matrix (II 1> . the corresponding kinetic parts of 
Eq. ( 12 II are given as 

T,| = -1-JdvJ dV'5(v-v')d 2 z a(v,v'), (22) 

Tl = -7^ JdV 5(r-r') Vi a(v, r'), (23) 
where = d 2 +d 2 , and the potential-energy contributions read 

V\\ = T ^f JdVz 2 a{v,v), (24) 

y± = rnuP_ j dv[x 2 + y 2 )a{v ^ (25) 

A change of the average energy due to the application of a single step of stochastic 
cooling is dominantly generated by the energy exchange with the externally applied optical 
fields that implement the momentum shift of atoms. The unidirectional flow of energy from 
the system to the optical fields is determined by the irreversibility introduced in the quantum 
measurement process. Apart from that, however, there are several sources of heating, among 
them also the back-action noise of the measurement itself. Especially at ultralow temperatures 
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a decay or even reversal of the net energy flow may be expected, when the detrimental heating 
terms compensate the sought cooling effect of the momentum shift. In the following we will 
extract all these energetic terms by considering the feedback-induced change of energy based 
on Eqs ( !2U -( l25t . More specifically, we consider the change of the average energy, i.e. the 
difference between the energy after a single step of stochastic cooling and that before, 

AE = E + -E_ = AT» + ATx + AVli + AV ± . (26) 



3.1. Energy change in the longitudinal motion 

The dominant change of energy will occur in the potential and kinetic energies associated 
with the longitudinal motion in z direction. The longitudinal potential energy Vji is given by 
Eq. 1241 and together with Eq. ( II 91 it can be shown that the change of longitudinal potential 
energy is 

Bid) 2 , - , 

= 8^" (27) 

This positive energy contribution arises from the back-action noise of the total-momentum 
measurement in z direction. The centre-of-mass z coordinate of the affected atoms, associated 
with the total mass m(N w ), is then subject to an increased uncertainty of the size (2(t) _1 , 
which Eq. J27I shows to introduce a heating term in the potential energy. 

For the kinetic energy of the motion in z direction a more involved calculation of 
the second-order z-derivative of the density matrix J19i is required. After some lengthy 
but straightforward calculation the following expression for the feedback-induced change of 
longitudinal kinetic energy is then obtained: 

AT„ = ^ <^_ + <^) 

2m(N w ) 2m(N w ) 2m(N w ) 2 

Here AN W = N w — (N w ) is the fluctuation of the actual atom number in the control beam 
around its average. The first term in Eq. d28b is the kinetic energy left in the system by the 
imprecise total-momentum measurement with resolution a. The second term is the sought 
cooling effect, where the centre-of-mass kinetic energy of the affected atoms is removed from 
the system. The last term, though, arises from quantum fluctuations of the number of atoms in 
the control-laser beam. This heating term appears since in the momentum-shift operation the 
actual atom number N w is not known but only estimated by (N w ). Thus this term represents 
a quantum-statistical imperfection of the feedback loop of stochastic cooling. 



3.2. Energy change in the transverse motion 

At first sight one may guess that the transverse motion is not affected by the feedback loop, 
since only momentum in z direction is measured and shifted. However, since the atoms 
that contribute to the measured signal are confined within the laser beam waist w±(r), the 
measurement of P w also contains an indirect measurement of the transverse position of atoms 
with a resolution roughly given by the diameter of the beam. Thus a back-action noise in 
the transverse momenta can be expected that may lead to further contributions to the kinetic 
energy. It is now left to show how large these energy contributions are compared with those 
emerging from the longitudinal motion. 

From Eqs dl9l and (125 \ it can be easily seen that the potential energy in transverse x and 
y directions is unchanged, i.e. 

AVx = 0. (29) 
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This result is obvious since only the momentum in z direction with a transverse spatial 
confinement is measured without affecting the noise in the transverse coordinates. Let us 
therefore consider the kinetic energy of the transverse coordinates as defined in Eq. d23t . 
Calculating the required second-order derivatives of Eq. ( 1 1 91 and performing the integrations, 
after some lengthy but straightforward calculus, we obtain for the change in transverse kinetic 
energy 

AT± = i/ ^[Vm(r)] 2 {^(fn(r)) + [ ff V + 3 w 2 ± (r)] tft (r) $ (r)) 

2cH (Av,)- 

-u, ± (r)({p z (r), P^}) 1 - - 1 „ - / dV z V Mi (r) ■ ({p(r), P w }), 



Aa 2 (N w ) " ) 2m(N u 

where {A, B}=AB+BA is the anti-commutator and the momentum density reads 

p(r) = - l - {^(r)V«£(r) - [vft(rj\ 0(r)} . (31) 
For a thermal equilibrium state, space dependent averages will have a symmetry with respect 



to z—> — z, and thus the second integral in Eq. d30i can be shown to vanish as an odd moment 
of z. 



3.3. Gaussian beam-waist function 

At this point the specific form of the beam-waist function shall be introduced. We consider 
here a Gaussian beam with the following definition for w± (r): 

(x-x ) 2 + (y~y ) 2 ~ 



w±(r) = exp 



In this way the area A^ of the integrated beam intensity, 



A 



dA w"]_ (r) = 7rr ( 



o- 



(32) 



(33) 



allows us to interpret ro as an effective radius of the laser beam. The squared gradient of 
w± (r) results from Eq. d32t as 



[V»i(r)f=<(r) 
and using these results, Eq. (I30i reads 



2 r„\ (x-x ) 2 + (y-y f 



r 4 



(34) 



1 



2a 2 <A^) 2 
1 



K^ 2 + >1W] ({^(r)0(r),P 2 }) 

(x-x a ) 2 + (y-y ) 2 



w ± (r){{p z (r),P w }} 



(35) 



This kinetic-energy change will determine the heating effect due to the transverse confinement 
of atoms in the laser beam. 



Transverse confinement in stochastic cooling of trapped atoms 



9 



4. Non-degenerate atomic vapour 

For ultracold temperatures near the condensation temperature Tq the longitudinal energy 
change has been discussed already in the approximation of a rectangular shape of the beam 
waist in Ref. [28]. In the following we evaluate the complete energy change J26l >. with the 
contributions given by Eqs. J27i . d28l >. i29i and J35i . We consider the regime of a non- 
degenerate gas, where the thermal de Broglie wavelength is much smaller than the interatomic 
distance. In that case the feature of indistinguishability of atoms can be neglected, keeping 
however the full wavemechanics of the single atom. That is, the single atom's position and 
momentum still obey the canonical commutator relation, from which several important effects 
emerge. 

Calculating expectation values for a thermal state in the canonical ensemble at 
temperature T and total atom number N, treating the atoms as distinguishable particles, the 
number of atoms in the laser beam, for example, results as 

<^>=^exp(-^). (36) 

Here we used the scaled distance d of the control beam from the trap origin and the scaled 
beam radius s, defined by 

d = \j { x l+vl)l L th, s = r /L th . (37) 
The rms extension of the atomic cloud is given by 

-1/2 

(38) 

with A.t = y/l/(2mu>) being the ground-state position uncertainty in the trap potential and 
fee being the Boltzmann constant. In the following we also use the size of the atomic cloud in 
units of the ground-state uncertainty: 

l th = L th /Ax . (39) 



L th = Ax 



tanh 



2fc B T 



4.1. Longitudinal energy change 



The complete change of energy in the longitudinal motion in units of vibrational energy quanta 
reads 



AE|/w 



(N w ) 



(<r/Ap ) 2 



/2 

AN 



(2 + s 2 ) 2 
s 2 (4+,s 2 ) 



exp 



/2 

_ Jk 
4 

Ad 2 



(40) 



(2 + S 2 )(4+s 2 ) 



- 1 



where ApoAxo = \ so that Apo = i/mw/2. The first term represents the measurement- 
induced noise leading to an increase of kinetic and potential energies. Taking into account 
only the energy change as given here this heating effect can be minimised by adapting the 
measurement resolution to the number of atoms in the control-laser beam as 



o- = Ap d(N, 



(41) 



The minimum heating due to this noise results then as one half energy quantum. 

The sought cooling effect is represented by the second term in Eq. (I40> . which is the 
centre-of-mass kinetic energy of the atoms addressed by the feedback. This value does neither 
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Figure 2. Boundaries between cooling (-) and heating (+) for the temperature T = 10 To and 
varying total numbers of atoms, cr is chosen as the optimal value given in Eq. EJ. 



depend on size nor location of the feedback region. For large temperatures Z t h - ► \J Zk^T / uj , 
so that the subtracted kinetic energy reduces to k^T /(2uj), manifesting the removed energy 
as being given by the equipartition theorem. 

Finally, the third term in Eq. J40I represents quantum noise due to the transverse spatial 
confinement of atoms subject to the feedback. It is produced by atom-number fluctuations 
that crucially depend via the average atom number on the size and location of the control- 
laser beam. This noise is always positive, leading thus to an unavoidable heating contribution, 
and does vanish only for s — -> oo. The limit s — > oo realizes the situation where all atoms are 
inside the control-laser beam, thus containing exactly N atoms with vanishing atom-number 
fluctuations. Moreover, the strength of this heating term diminishes with increasing total 
number of atoms, as can be observed in Fig. [2] 

There the contours A£7i| = have been plotted in the parameter space (s, d) for varying 
total atom numbers at fixed temperature T= 10 To, where 



is the condensation temperature in the thermodynamic limit [40] with ((n) being the Riemann 
£ function. These contours represent the boundaries between cooling (-) on the right-hand side 
of the contour and heating (+) on its left-hand side. They show that finer spatial resolutions 
s and larger distances from the trap centre d are allowed for cooling when the total atom 
number increases. Whereas the measurement-induced noise is constant and negligibly small, 
the major contribution leading to a restriction of the parameters s and d comes from the atom- 
number fluctuations in the chosen beam waist. 




(42) 
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Figure 3. Contours A_Ey 
temperatures and N = 10 6 . 



(dashed lines) and AE = (solid lines) for different 



4.2. Transverse energy change 

The energy change in the transverse motion can be obtained in the same way as that for the 
longitudinal one, cf. Eq. J40i . and reads 



AE ± /u 



N 



(N w ) , (a/A P o) 2 



(a/Ap ) 2 (N w ) 



N 



4{N U 



i 2 +r 2 

'th ^ 'th 



(N w ) (v/ApoY 
I 2 

'th 



N(N-l)l 2 h s 2 [A + s 2 (2 + d 2 )] 



4 + s 2 (2 + rf 2 ) 
(2 + s 2 ) 3 

+ s 2 (2 + ri 2 ) 
(4 + s 2 ) 3 

2d 2 



cxp 



■ exp 



2 + s 2 
2d 2 



4 + s 2 



( Nw ) 2 4 (2 + S 2)4 
1 



cxp 



2 + s 2 



4( ( x/Ap ) 2 (7V u 



12 + s 2 ( 2 + d 2 ) 
— — — - exp 



(6 + s 2 ) 3 



3d 2 



1 \ s 2 8 + s 2 2 + d 2 

+ iV(iV-l)- ? — v 9 .„ ; exp 

v '2 + s 2 (4 + s 2 ) 3 



6 + s 2 

d 2 (8 + 3s 2 ) 
(4 + s 2 )(2 + s 2 ) 



(43) 



As mentioned before these heating contributions vanish in the limit s — > oo. Moreover, they 
depend on the total number of atoms N, which is also mediated by the dependence on (N w ), 
see Eq. (I36> . and possibly on a. 

Let us first consider the emerging changes in the boundary between cooling and heating, 
when now also the transverse energy change is taken into account. That is, we look for the 
contour in (s, d) parameter space that satisfies AE = 0, where AE = AE\\ + AE±. In Fig. [3] 
this contour (solid curve) is shown and compared to the boundary based only on longitudinal 
terms (dashed curve) for a total atom number of N = 10 6 . It is clearly observed that the 
additional heating terms due to the transverse confinement of atoms in the control beam leads 
to a shift of the boundary to larger values of s and smaller values of d. The former is due 
to the fact that the measurement of the total momentum in z direction of atoms in the beam, 
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Figure 4. Boundaries AE = for T = 10 To and varying total atom numbers: TV = 1 (dotted 
curve), N = 10 (dashed curve), N = 100 (dot-dashed curve), N = 10 6 (solid curve). 



indirectly also represents a measurement of transverse coordinates within the beam-waist size. 
This leads again to measurement back-action noise, which is naturally reduced by increasing 
s. 

The fact that the boundary is shifted to smaller values of d is due to relative atom-number 
fluctuations that decrease for increasing atom numbers found near to the trap centre. It should 
be noted that either curve actually contains two different curves at temperatures T = 10 T 
and T = 10000 T , which however cannot be distinguished in the plot. The only explicit 
dependence in Eqs. d40l and d43t is due to the occurrence of l t ^. For JV 3> 1, however, 
this dependence is very weak. Nevertheless, all features discussed here and in the following 
implicitly depend on temperature via the chosen scaling of s and d by L t h- 

Assuming that for increasing total atom number the optimised measurement resolution 
g increases, as for the value given in Eq. (14-li for longitudinal terms, for large atom numbers 
N 3> 1 an asymptotic expansion can be given for Eq. (I43> : 



AEj_/uj = - 
4 



(N w ) | (<x/Apo) s 
(a/Apo) 2 (N w ) 



+ l 2 



th 



4+s 2 (2+d 2 ) 



s 2 (2 + s 2 ) 2 



(44) 



where Eq. d36l > has been used. In the same limit the longitudinal energy change can be 
approximated to finally obtain the total change of energy as 



AE/uj 



(N w ) 



{a/Apo 



1- 



4+s 2 (2 + rf 2 ) 



s 2 (2 + s 2 ) 2 



/2 

Hh 



1- 



4+s 2 (2+rf 2 ) 
s 2 (2 + s 2 ) 2 



■(45) 



This asymptotic expansion represents a good approximation starting already from atom 
numbers N > 100, as can be seen from Fig. [4] There it is observed that for fixed temperature 
the boundaries converge quickly to the asymptotic one for N — > oo. 

Equation i45\ shows that again the optimal value for the measurement resolution is given 
by (14 1 i . for which the final result reads 



AE/u 



4+s 2 (2 + rf 2 )' 
s 2 (2 + s 2 ) 2 _ 



/2 
'th 



1 - 



4 + s 2 (2 + d 2 )' 
s 2 (2 + s 2 ) 2 _ 



(46) 
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Figure 5. Dependence of s m m 011 l?h f° r large total number of atoms N — > oo. 



For this expression the solution of the condition AE = can be analytically given as 



d(s) = (2 + s 2 ) 1 



I 2 + 2 



s 2 (2 + s 2 ) 



(47) 




For d = the corresponding value for s is the minimal beam-waist radius. This minimal radius 
is obtained as 



(48) 



In Fig. 13 this function is shown in dependence on l 2 h , which for high temperatures is 
proportional to T. At l 2 h = 2 the removed centre-of-mass kinetic energy exactly compensates 
the measurement-induced heating, which requires s m j n — ► oo. For values l 2 h < 2 cooling does 
not occur, since the unavoidable measurement-induced noise can no longer be compensated. 
However, the corresponding temperatures for l 2 h < 2 are below the condensation temperature 
To, where our approach for a non-degenerate gas is no longer valid. 
Ref. [28]. 



For this regime see 



For larger values of l 2 h , or T correspondingly, the limiting value s m [ n — * (vo— l) 1 / 2 m 
0.86 is reached. In this regime only a fraction of the atomic cloud needs to be subject to the 
feedback loop since s m j n < 1, The unsealed minimum beam-waist radius is thus 86% of the 
rms extension of the atomic cloud. 



5. Conclusions 



In summary we have studied the effects of transverse confinement in stochastic cooling of 
trapped atoms. It could be clearly shown that these effects are substantial for the cooling 
process and that minimum values for both the size and location of the control-laser beam 
exist. 

In the regime of non-degenerated gases analytical expressions could be derived, 
that contain the full quantum-fluctuation effects. Among these effects are atom-number 
fluctuations that appear due to the finite volume of the control beam. They appear in form 
of an imperfection of the feedback loop and in form of back-action noise due to the indirect 
measurement of transverse coordinates. 
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